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Abstract 

Within the context of supersymmetric hybrid inflation, we study a multiple-stage infla- 
tionary scenario that can generate a primordial spectrum of adiabatic density perturbations 
with a break at kh ~ 0.05 h Mpc~^. The presence of such a break is supported by the APM 
galaxy survey data. We consider a specific model within this scenario and confront it with 
observational data. We reproduce the angular power spectrum of CMB anisotropics, and 
also account for the break in the power spectrum of galaxy clustering. In addition, we find a 
value for cts in agreement with the one deduced from observations. A characteristic property 
of the spectrum is a drop of the spectral index from n ~ 1 to n ~ 0.6 at fc ~ fcf,. 



PACS number: 98.80.Cq, 98.80.-k, 98.95.D 



1 Introduction 



The origin of the measured anisotropies in the cosmic microwave background and the gener- 
ation and evolution of large-scale structure in the universe, such as the observed clustering 
in the galaxy distribution, represent outstanding questions in modern cosmology. Within the 
framework of gravitational instability, there are two currently investigated families of models 
to explain the formation of the observed structure. Initial density perturbations can either be 
due to "freezing in" of quantum fluctuations of a scalar field during an inflationary period [j^], 
or they may be seeded by topological defects, which can form naturally during a symmetry 
breaking phase transition in the early universe Q. The cosmic microwave background (CMB) 
anisotropies provide a link between theoretical predictions and observational data, which may 
allow us to distinguish between inflationary models and defect scenarios by purely linear analysis 

i- 

Within the inflationary paradigm, a possible mechanism for the generation of the angular 
power spectrum of CMB anisotropies and the creation of the large-scale structure is based on 
the quantum fluctuations that exited the horizon during inflation Q . These fluctuations are the 
source of the primordial spectrum of density inhomogeneities |^, which has left an imprint on 
the CMB. The observed large-scale structure could have been generated by the growth through 
gravitational instability of this primordial spectrum of perturbations in the otherwise uniform 
distribution of matter. 

By making some assumptions about the matter content of the early universe, one can in- 
vestigate whether a given spectrum of primordial perturbations can evolve into the large-scale 
structure observed today. The evolution of the spectrum depends crucially on the nature and 
amount of the dark matter in the universe. The two extreme cases, of hot dark matter (HDM) 
and of cold dark matter (CDM), lead to inconsistencies with observational data. The HDM 
model has been abandoned, since the thermal motion of massive neutrinos wipes out any small- 
scale structure On the other hand, the standard CDM model with a flat (spectral index 
n = 1) initial spectrum of adiabatic perturbations predicts too much power on small scales, if 
one normalizes it to the COBE data at large scales. The standard CDM model leads to a value 
(Tg ~ 1.2 for the variance of the total mass fluctuation in a sphere of radius 8 Mpc//i Q. On the 
other hand, the observational value, inferred from the abundances of rich clusters of galaxies, is 
as -0.6 ±0.2 |, |. 

As the standard CDM model is a generally successful model (e.g., it successfully reproduces 
galaxy clustering statistics, the various epochs of structure formation, as well as peculiar velocity 
flows, while it maintains CMB anisotropies at a level lower than observational upper limits), 
several approaches have been proposed, which modify its underlying assumptions (for a review 
see Ref. In order to render CDM models compatible with observational data, one possibility 

is to change the initial spectrum of perturbations by introducing a tilted (spectral index n < 1) 
scale- free spectrum. However, this approach has not led to successful predictions [0]. Another 
approach is to consider less standard values for the cosmological parameters, while leaving the 
initial spectral index unchanged (i.e., n ~ 1). For instance, one can either investigate CDM 
models with a higher fraction of baryons ll^ , or models which have a mixture of CDM with an 



amount of HDM, or even models with a positive cosmological constant 14]. 

Another interesting possibility is to consider power spectra with a break, which arise natu- 
rally in models with more than one stages of inflation - an approach adopted by several authors 



|5|]-|18]. A model of double chaotic inflation which generates such a spectrum has been studied 



in detail in Refs. |]1^, 19, 20 1. The model contains two scalar fields with disconnected quadratic 
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potentials. 

In this paper, we discuss the possibiUty of a multiple-stage inflationary scenario within 
the context of supersymmetric hybrid inflation [^]. The necessary superpotential can arise in 
supersymmetric GUT models. Moreover, the observable part of inflation takes place for field 
values below the Planck scale, so that supergravity corrections are under control. Inflationary 
models such as the one we are considering here have been discussed in Ref. and have been 
shown to survive the supergravity corrections in Ref. |23]. 

A strong motivation for these models has its origin in the issue of the initial conditions 
- p^ that are necessary for hybrid inflation [^] |^. This is a consequence of the presence 
of (one or more) scalar fields orthogonal to the inflaton, and the constraint from the COBE 
data on the inflationary energy scale V^^'^ (determined by the vacuum energy density during 
inflation), which must be at least two or three orders of magnitude smaller than the Planck 
scale 0. The difference between the energy scale mpi at the end of the Planck era, when 
classical general relativity becomes applicable, and the inflationary scale V^^^ implies that the 
various fields evolve for a long time before settling down along the inflationary trajectory. The 
Hubble parameter H sets the scale for the "friction" term in the field evolution equations. 
When the energy density drops much below rripi, the smallness of this "friction" term results 
in a very long evolution, during which the fields oscillate around zero many times. Some of 
the trajectories eventually settle down in the valley of the potential that produces inflation. 
However, the sensitivity to the initial conditions is high because of the long evolution. A slight 
variation of the initial field values separates inflationary trajectories from trajectories that lead 
to the minima of the potential, where inflation does not occur. As a result, the initial field 
configuration must be extremely fine-tuned for inflation to set in. The fields orthogonal to the 
inflaton (whose typical fluctuations at the end of the Planck era are ~ mp;) must be zero with 
an accuracy of at least ~ 10~^mpi over regions that exceed the initial Hubble length by one or 
two orders of magnitude . 

In Ref. a simple resolution of the issue of flne-tuning described above was suggested. A 
scenario with two stages of inflation was proposed within the context of global super symmetry. 
The first stage has a typical scale ~ mp;, which implies that the "friction" term in the field 
evolution equations of the fields is large. As a result, the system settles down quickly along an 
almost fiat direction of the potential and this stage of infiation occurs naturally. By generating 
an exponential expansion of the initial region of space, it also provides the homogeneity that 
is necessary for the second stage. The latter has a characteristic scale much below mpi and 
generates the density perturbations that result in the CMB anisotropy observed by COBE. This 
scenario was generalized in the context of supergravity in Ref. |2^] . 

It is a logical next step to contemplate the possibility of a multiple-stage inflationary scenario, 
with a sequence of scales that starts near mpi and continues down to the scale implied by COBE, 
or even to lower energy scales. For inflation driven by one fleld in the context of supergravity, 
such a possibility has already been proposed in Ref. [18| and has been associated with a break 
at kb — 0.05 h Mpc~^ in the power spectrum of galaxy clustering derived from the APM galaxy 
survey. This break has been linked to a possible phase transition between two of the inflationary 
stages. In this paper we would like to provide another realization of such a scenario, within a 
multi-field theory, with an estimation of the predicted power spectrum. The determination of 
the spectrum during the phase transition is beyond our technical capabilities at the moment. 
However, the phase transition affects only a small range of scales, while the rest of the power 



^An alternative mechanism for the resolution of this problem has been proposed in Ref. Pq ]. 
^Throughout the paper we use the "reduced" Planck scale mpi = Mpj/VStt, Mpi = 1.22 x 10^^ GeV. 
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spectrum is strongly constrained by the observational data. 

Firstly, we are interested in examining whether the COBE data can be made compatible 
with the data of galaxy surveys, employing the smallest number of cosmological parameters. 
More specifically, we consider a CDM model, with zero cosmological constant, present Hubble 
parameter h = 0.5 {Hq = h 100 km/s/Mpc) and flat geometry ^matter = 1- We assume that 
the fraction of the critical density in baryons is Of, = 0.05 and in cold dark matter ^cdm = 
0.95. Subsequently, we repeat the calculation, allowing for a non-zero value of the cosmological 
constant, a scenario supported by recent observations |14|. 

In Section 2 we discuss the general form of the power spectrum predicted by the APM data. 
In Section 3 we present our model and in Section 4 the properties of the inflationary stages. 
The primordial power spectrum in our model is calculated in Section 5 and is compared with 
observations in Section 6. Our conclusions are given in Section 7. 



2 The APM power spectrum 



An indication for the possibility of more than one stages of inflation is provided by the power 



spectrum deduced from the APM galaxy survey |30|. This spectrum has a break at a charac- 

-1 



One explanation for this break could be a 



teristic scale fc^ ~ 0.05 h Mpc [31 



18, 



sharp change of the spectral index of the primordial spectrum between two stages of inflation 
|18|. In this section we summarize the main points of the analysis in the above references, in 
order to motivate the multiple-stage scenario. In the following sections we shall show how such 
a scenario can be realized in the context of a specific model. 

In Fig. 1 we plot the logarithmic slope of the power spectrum recovered from the APM 
angular galaxy catalogue [31, (We have used the data of Table 2 in Ref. [|3S 



33[. (We have used the data of Table 2 in Ref. |33||.) The slope 
has been calculated through a least-squares fit for groups of 3 data points, assuming that the 
spectrum is locally linear in a logarithmic scale. The results are in agreement with the detailed 
analysis of Ref. [^], even though our estimated errors are slightly different. A sharp change in 
the slope of the spectrum is visible around /c^ ~ 0.05 h Mpc^^. 

This break persists even after the effects of the non-linear evolution of matter fluctuations 
(when the density contrast becomes of order 1) have been removed. Through A^-body sim- 
ulations, Baugh and Gaztahaga [Esl have estimated the linear spectrum that gives rise to the 



observed spectrum for a spatially flat universe with critical matter density and zero cosmological 



constant. Their estimate is well fitted by the curve [33, 34] 



(2.1) 



with C ~ 7.0 X 10^, kb ~ 0.05 h Mpc^^, a ~ 1, /3 ~ 1.6. This result is in agreement with 
the spectrum obtained through the application of the formula of Jain et al. [ p^] f or the relation 
between the linear and non-linear spectra. The formula of Peacock and Dodds |3g] does not lead 
to very good agreement with the results of the A^-body simulatio ns p^ ]. In Fig. 1 we display 
the logarithmic slope of the linear spectrum as predicted by Eq. ( |2.lD (solid line). The sharp 
decrease in the slope is again apparent. 

The implications for the primordial spectrum can be obtained if one deconvolutes the spec- 
trum by using the linear transfer function, which determines the scale-dependent growth of 
linear perturbations. A parametrization of this function is given by the relation 



T{k) = [l + [ak + (6A:)3/2 ^ (^^^2^ 



-l/v 



(2.2) 
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with a = 6.4 T-^h-^Mpc, 6 = 3 T-^h-^Mpc, c = 1.7 r^^^-^Mpc, u = 1.13 and T = Qhe''^^'', 
where Of, is the fraction of critical density in baryonic matter. For h = 0.5, = 0.05, the 
spectral index n of the primordial spectrum P{k) ~ Pi{k)/T'^{k) is given by the dashed line in 
Fig. 1. A transition from a value n ~ 1 at large scales to an average value n ~ 0.5 at small scales 
is observed. It must be pointed out that the above analysis cannot be considered conclusive, 
as the parametrization of Eq. ( |2.2| ) may not be sufficiently accurate. In section 6 we present a 
detailed comparison with the linear power spectrum of Eq. ( |2.1| ) using the code CMBFAST of 
Seljak and Zaldarriaga [^ ]. 

The variation of the spectral index can be explained by the generation of density perturba- 
tions during two stages of inflation, driven by fields with different potentials. The first stage 
must have lasted long enough for scales between k ~ 3 X 10"^ h Mpc~^ (relevant for the COBE 
data) and /c ~ /c^ ~ 5 x 10~^ h Mpc~^ to have crossed outside the horizon during its dura- 
tion. This requires at least 5 e-foldings to be generated by the first stage of inflation. The 
large deviation of the spectral index from 1 has strong implications for the duration of the 
second stage. For the models of interest to us, the inflationary potential has the general form 
V ~ Vo(l -|- cf{(f)/mpi)), where c is a small constant. Moreover, the inflaton fleld is constrained 
to he (j) ^ mpi. A general analysis of the spectral index predicted by such potentials, for various 
choices of the function f{(j)/mpi), is given in Ref. [^^. According to Table 1 of that reference, 
the total number of e-foldings is constrained by the relation |n — l|(A^/50) ^ 0.08. For n ~ 0.6 
we infer N ^ 8. More specifically, for the case f{(j)/mpi) = \n{(f)/mpi) that is relevant for our 
model, |n - l|(A/50) ~ 0.02 and AT ~ 2. 

The above considerations support a picture of multiple short bursts of inflation, with signifl- 
cant variations of the spectral index of the primordial spectrum. The observational data provide 
information on two of these inflationary stages. The COBE-DMR measurements require a stage 
with n ~ 1 and at least 5 e-foldings, while the APM galaxy survey data support the possibility 
of a second stage with n ~ 0.6 that generates ~ 2 e-foldings. The total number of e-foldings 
during all the inflationary stages with observable consequences must be ^ 60 for the flatness 
and horizon problems of standard cosmology to be resolved. However, the stages responsible 
for the last ~ 50 e-foldings generate density perturbations at scales that are strongly affected 
by the non-linear evolution. As a result, it is very difficult to extract information on their 
characteristics. 

The above picture contrasts with the double inflationary scenario of Kefs, m, 1^, where 
inflation is driven by two disconnected quadratic potentials for the flelds ui, a2- The first 
stage of inflation occurs for ai^2 ^ "ipi and terminates when ai ~ mpi, while the second one 
terminates when a2 ~ mpi. For the second stage to generate less than ~ 60 e-foldings, so 
that the break in the spectrum is observable, the value of (T2 at its beginning must not be very 
far from nipi. In our model there is no need for a similar adjustment. The requirement of a 
spectral index signiflcantly below 1 for the part of the spectrum generated by the second stage 
of inflation automatically constrains the number of e-foldings from this stage to be less than 60, 
independently of the value of a2- 

In the following section we describe a model which provides the potential for the two stages 
of inflation that are relevant for the COBE and APM data. 
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3 The model 



The model we consider is described by the superpotential 

W = Si (-^? + \i^i<^i + g^2^2) + S2 + A2^2$2) ■ (3.1) 

The superfields $1, $1 and $2) ^2 transform under internal gauge symmetry groups A4i and 
M.2 respectively. The Si, S2 superfields are gauge singlets. We have assumed that ^2 is the 
linear combination of gauge singlets that does not couple to |2^. We make the simplifying 

assumption that the scalar components of the various superfields are real. The imaginary com- 
ponents do not alter the qualitative picture of inflation in this model, while they would make 
the calculation of the spectrum of density perturbations much more complicated. Staying along 
the D-Qat directions, we define canonically normalized scalar fields according to 

^2 ^2 = ^2= ^. (3.2) 

The potential is then given by the expression 



V{ai,a2,4>i,4>2) = 



dW 



2 i2\^ 9 2i2 , 9^ lA , 9^1 ,2 i2 , "^1 2i,2 

+ (f^l-^4) +7 (5^1 + ^2^x2)' 02, (3.3) 



where the mass scales /ii, /i2 are chosen fii > ^2- The minima of this potential are located at 

n i2 ^ 2 ^9 2 i2 ^ 2 l<i A\ 

0-1 = 0-2 = 0, 01 = — /ii - T-^^2, 92 = ^/^2- (3-4) 

Ai A1A2 A2 

For = (/)2 = the potential is independent of 0-1,2- Its value V = fif + ^2 gives the vacuum 
energy density during the first stage of inflation. Along this direction, the mass terms of the 
(pi^2 flelds are 

Ml = - Ai/x? + (3.5) 
Ml = - X2f4 -gA + \ (9(^1 + ^2(72? ■ (3.6) 
The mass term Ml becomes negative for 

Al 

This indicates the presence of an instability which can lead to the growth of the (pi field. 



We use the same notation for the superfields and their scalar components. 
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Another flat direction (independent of (T2) exists for ai = 0, = 4/i^/Ai, = 0. The value 
of the potential V = ^2 gives the vacuum energy density during the second stage of inflation. 
The mass term of the 02 field is given by 



Mi 



\ ,,2 I ^^2^2 
^2^2 + Y^2- 



An instability appears for 



2 ^ 2 

'^2 ^ '^2ins 



2^4 

\2 ' 



(3. 



(3.^ 



which can lead to the growth of the 02 held. 

The flatness of the potential is lifted by radiative corrections. During the first stage of 
inflation and for ai^2 far above the instability points, the one-loop contribution to the effective 
potential is ^ 



Ay(cri, 0-2) ^ 



Ml 
167r2 



Allll 



2A2 



In 



+ 



+ 



M2 
167r2 



{>^2l4 + 5^1 



In 



(ffUl + A20-2) 

2Ai 



+ 



(3.10) 



where Mi^2 are the dimensionalities of the representations of the groups A^i,2 to which the 
superfields ^1^2 belong. The exact value of the normalization scales A1.2 is not important for 
our discussion. In the following we shall consider theories in which the two sectors (5i,$i) 
and (5*2, $2) are essentially decoupled. For this reason we shall assume that the coupling g is 
sufficiently small for the radiative correction of Eq. ( |3.1C1| ) to be approximated as 



In 



AM 



M2 f. 2 2\2 

+ 16^i^2''2 + 5ftJ 




(3.11) 



The above contribution provides the slope that leads to the slow rolling of the ai^2 fields during 
the first stage of inflation. 

During the second stage of inflation the slope for the (T2 field is provided by the radiative 
correction 



Ay(c7i,a2) 



M2 
167r2 



A2/^2 



In 



2A?, 



+ 



(3.12) 



Away from the flat directions the radiative corrections are small and we neglect them. 

An important question concerns the supergravity corrections to the above potential. In par- 
ticular, it is important that the flat directions of the potential are not lifted by these corrections 



|21|. A detailed analysis of this problem is presented in Ref. |23], where it is shown that, for 
the flat directions to be preserved, at least one of the two inflationary stages must be driven 



by D-term energy density |39|. However, the form of the potential is the same as the one we 
described above. Two flat directions exist with a small slope generated by logarithmic radiative 
corrections. 



It must be pointed out that the supersymmetric cancellations that lead to Eqs. (3.10), (3.12) require the 
inclusion of radiative corrections from both the real and imaginary scalar components of the superfields "l>i, $2. 
However, this is the only effe ct of the imaginary components on the infiationary stages. This is the reason why 
we omitted them in Eq. (3^). For a discussion that includes these components see Ref. |22| . 
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4 The inflationary stages 



The model we described in the previous section allows for two stages of inflation separated by 
an intermediate stage. 



4.1 The first stage of inflation 

The Hubble parameter during the first stage of inflation is almost constant 

Hi ^ (4-1) 

We assume that the part of this stage with observable consequences starts at a time tu and 
finishes at a time tij. Because of our assumption > fi2, the vacuum energy density associated 
with the fields ai, 0i dominates the first stage. Therefore, this stage terminates when the "slow- 
roll" conditions for cji are not satisfied any more. For the potential of Eqs. ( |3l3|) , ( |3.11|) , this 
happens for 

^!{tif)^CTjf^^mli. (4.2) 

Unless the coupling Ai is taken much smaller than 1, cr^j is much larger than the instability 
point aji^, of Eq. (|3?7|) . 

The number of e-foldings from a time t until the end of the first stage is given by 

\a{t) J MiXf /if rripi 

where a{t) is the scale factor and ai/ = a(ii/) • The total number of e-foldings during the 
observable part of this stage is given by 

„,„=,„(M)=j!]l*l±f4f?l^, (4.4) 

\auJ MiXf /xf mf,i 

with cr^j = af {tu), au = a (tu). The relative change of the values of the two fields during the 
first stage is 

In the following, we shall use parameters that make this ratio much smaller than 1, so that the 
evolution of o"2 can be neglected during the first stage of inflation. 



4.2 The intermediate stage 

The intermediate stage lasts between times tif and t2i- After the end of the slow-roll regime, 
the o"! held quickly rolls beyond the instability point of Eq. ( |3.7D . Subsequently, large domains 
start appearing in which the value of the (pi field grows exponentially with time. For statistical 
systems, for which the expansion of the universe is not relevant, this process is characterized as 
spinodal decomposition. The expansion of the universe complicates the above picture, but the 
details are not important for our discussion. We assume that this initial stage of instability is 
fast, and soon the fields take values away from the ai axis and in the vicinity of the minimum 
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at cJi = 0, (pi = 4^1/ Xi, where the curvature of the potential is positive. Our assumption is 

reasonable because the ai field rolls to the origin within a time ~ H^^ ~ \/3mpi/^ fif + //^ or 
slightly larger. On the other hand, the typical time scale for the growth of the (pi field is given 
by the absolute value of the curvature at the origin and is ~ (\/AiA*i) • As a result, we expect 
that (pi grows to a value near the minimum within a fraction of a Hubble time. 

After a short complicated evolution, the massive fields ai, (f>i settle into a regular oscillatory 
pattern around the minimum, with the universe characterized by an equation of state p = wp. 
For a system of massive oscillating fields, such as the one we are considering, zz; = 0. Prom this 
point on, the energy density of the oscillating fields is dissipated through the expansion of the 
universe. If the fields have decay channels into lighter species that eventually thermalize, the 
equation of state of the radiation-dominated universe has w = 1/3. When the energy density 
becomes comparable to the second stage of inflation can begin. During the intermediate 
stage, the scale factor increases by a total amount 

N,^, = lJ^\= In f ^ V (4.6) 

where a2i = a {t2i). 

The effect of the oscillations of ai, cpi on the stability of (T2, (/>2 is minimized if the former 
have fast decay channels into lighter species. Under this assumption, the fields (T2, (1)2 remain 
constant during the intermediate stage, even when the supergravity corrections are taken into 
account [^]. In the following we make this assumption and use u; = 1/3 during the whole 
intermediate stage. 



4.3 The second stage of inflation 

This stage starts at a time t2i and finishes at a time t2/- The Hubble parameter is 

Infiation stops when the "slow-roll" conditions for (T2 are not satisfied any more. For the potential 
of Eqs. (U), (I3.12D , this happens for 



M2A2 



a'2 {t2f) = ^ -^^Pi- (4-8) 

Unless the coupling A2 is taken much smaller than 1, fj^j is much larger than the instability 
point ofEq. (Q. 

The number of e-foldings from a time t until the end of the second stage is given by 

( a2f\ _ 4tt^ ^K^)-^!/ 



\a{t)J M2XI rn?pi 
where 02/ = o (^2/) • The total number of e-foldings during this stage is given by 

N2tot = In = ,0 o ^ (4-10) 



a2: 



J M2XI 



with a2i = cr| {t2i) 
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In the next sections we shall consider parameters for which the above equation predicts 
N2tot — 2 e-foldings. As the accuracy of the analytical expressions presented in this section is 
questionable for so small N2tot, we have verified the properties of the second inflationary stage 
numerically. We have found a value for N2tot that is close to 3 for our choice of parameters. 



4.4 The remaining stages of inflation 

The total number of e-foldings during all the inflationary stages with observable consequences 
must be ~ 60 for the flatness and horizon problems of standard cosmology to be resolved. This 
means that additional inflationary stages must take place after the first two that we discussed 
above. 

A simple extension of our model that would provide this possibility involves an additional 
term 5*3 (— /^l + As^a^I^s) in the superpotential. We assume that the new system of superfields 
5*3, ^3, $3 is decoupled from the ones driving the first two stages of inflation. The potential and 
the properties of the third inflationary stage are completely analogous to the ones we described 
above. For a sufficiently small scale fi^ the discussion of the previous subsections is not affected. 
A third stage of inflation starts after a second intermediate stage, as soon as the energy density 
becomes comparable to /if. The total number of e-foldings during the third stage is 

Nstot = In [ ^ )= ——^ 2-^' 4.11 

V a3i J M3AI mf,i 

where cr'^f/'rn'pi ~ M^sA^/Stt^ and the various quantities are defined in analogy to the first two 
stages. One can envisage several stages driven by such systems of superfields. The calculation of 
the supergravity corrections to the potential of such complicated systems of fields is a difficult 
task. However, we assume that it is possible to guarantee the presence of the necessary number 
of flat directions, even when supergravity corrrections are taken into account. 

There are other possible sources of inflation as well. The increase of temperature during the 
reheating after the first inflationary stages may lead to the restoration of spontaneously broken 
symmetries. Additional stages of inflation may take place when the temperature falls below the 
energy scale associated with these symmetries. Inflationary stages, such as the thermal inflation 



of Ref. |4C], may occur even as low as at the TeV scale. 



5 The primordial spectrum 

The primordial spectrum of density inhomogeneities has its origin in the quantum fluctuations 
that crossed outside the horizon during inflation |^] . The scale k crosses the Hubble radius when 
k = aH. The scale kcoBE/o-o ~ Hq exited the horizon at the beginning of the observable part 
of the first stage of inflation. This implies that uoHq = auHi. (Here oq, Hq ~ (3000 Mpc)~^ h 
are the present values of the scale factor and the Hubble parameter. Following the standard 
convention, we set ao = 1 in the following.) 

For k ^ k2 = a2iH2 the Hubble radius is crossed only once, during the first stage of inflation, 
while for k ^ ki = aifHi it is crossed once during the second stage. The scales k2 < k < ki, 
however, cross the Hubble radius three times: They exit the horizon during the first stage, re- 
enter duing the intermediate stage and exit again during the second stage. We define the scale 
K according to 

K=y/hh (5.1) 
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and we obtain 



-1 



1 1 / H 

Nitot + ^Nint + - In (^-^ ) , 



with Nitot, Nint given by Eqs. (4^), ( [4. 61) respectively. For the ratio k2/ki we obtain 

'A;2\ 2 



In 



(5.2) 



(5.3) 



The gravitational potential for scalar metric perturbations that are generated during the 



first stage of inflation is given by the expression [16, 41, 42| 



$= 2C— + D- 



m H2mli{Vi + V2y 



(5.4) 



with 



C 
D 



- \h ( 

' 2 [ 1^1 + ^2 CTi 

1 / 5ai 6a2 
3H \ &i &2 



+ 



V2 5(T2 
Vi + V2 CT2 



k=aH 



(5.5) 
(5.6) 



k=aH 



The first term in the right-hand side of Eq. (|5.4| ) can be interpreted as the adiabatic contribution 
to the spectrum, while the second one corresponds to entropic fluctuations. Until the end of the 
first stage of inflation (when ti — ti/ according to Eq. ( [4. 2D ), the functions Vi((Ti), V2(<T2) are 
the potentials of the two disconnected fields cJi, 02 and can be derived from Eqs. ( |3.3| ), (|3.11[ ): 



V2{a2) 



4 

^J'2 



+ 



167r^ 



M2 
167r2 



Afof 
2A2 



In 



+ 



(A2M2 + 9^^l 




(5.7) 



However, the subsequent evolution of the entropic contribution is very complicated, due to the 
growth of the (j)i field and the possible decay of o"i, </>i into lighter species. For this reason, 
we concentrate on the adiabatic contribution, which can be followed until today. We point 
out, however, that for our choice of parameters of the model, the total spectrum of adiabatic 
perturbations is dominated by the fluctuations of the (J2 field, with the contribution of <ti 
being negligible (see below). Since the entropic contribution arises because of the simultaneous 
fluctuations of at least two fields, the dominance of the (T2 fiuctuations indicates that the entropic 
contribution to the spectrum might be subleading. 

Assuming that the field fluctuations are random gaussian variables, we obtain for the spec- 
trum of adiabatic density perturbations 



5l{k) 



1 



257r^ 




(5. 



k=aH 



where primes on Vi(ai), ^2(0-2) denote derivatives with respect to cJi, C72 respectively. We 
assume that the evolution of o"2 is negligible during the first stage of inflation (so that 02 — 02% 
during this whole stage). This is guaranteed if the ratio of Eq. (4^) is much smaller than 1. 
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We also concentrate on the case Vi/V{ <C V2/V2, for which the metric perturbations are mainly 
generated by the fluctuations of the o"2 field during both stages of inflation. This requires 



Ml Af 



M2 \ \2i4 + gut 



and leads to 



Stt 



Hia2i 



5M2 \X2fii + gfii 



ni 



(5.9) 



(5.10) 



PI 



19, Eg], which has a logarithmic 



The predicted spectrum is scale invariant with a spectral index n = 1 to a very good accuracy. 
It provides a good approximation to the primordial spectrum for k < k2. It should be con- 
trasted with the spectrum predicted by the model of Refs. | 
dependence oc In^^'^ (K/k) for small k. 

During the second stage of inflation the metric perturbations are generated by the a2 held 
and the spectrum is 

Svr H2 [a2]k=aH2 



SH{k) 



5M2XI 



m 



(5.11) 



PI 



The fleld cj2 evolves from a2i to cj2/ during this stage and the resulting spectrum has a spectral 
index different from 1. The scale dependence is given by the relation 



0"2 



1 2 



mpi 



'2i 



+ 



k=aH2 



m 



PI 



M2XI 



Nuot + N, 



int 



In 



In 



Hi 
H2 



(5.12) 



The resulting spectral index is 



n 



1 



d\nk 



M2XI 



m 



PI 



4^2 a^(^ky 



(5.13) 



with cr2{k) given by Eq. (5.12). We are interested in scales that crossed outside the horizon 
during the beginning of the second stage of inflation. For them the spectral index is given by 
the above expression with (T2{k) ~ (T2i- 

The calculation of the form of the spectrum for scales k2 < k < ki that re-enter the horizon 
during the intermediate stage is prohibited by the extremely complicated nature of this stage. 
This is due to the presence of two massive oscillating and decaying flelds (ai^cpi). However, 



we do not expect any strong features (such as spikes |43|) in the spectrum associated with 
the intermediate stage. The situation is very similar to the "double inflation with a break" 
of Ref. [^]. In that case there is only one massive oscillating field with no decay channels, 
whose energy is dissipated through expansion. The detailed calculation of the spectrum has not 
revealed any strong features. In our model, no "massless" fields exist during the infiationary 
stages, other than cji, cj2 whose contribution to the spectrum we have computed. The first stage 
of inflation ends for cii ~ cJi/ » cruns, long before the orthogonal fleld (pi becomes massless or 
develops an instability; and similarly for the second stage. A smooth interpolation between the 
two parts of the spectrum that we have computed is likely, such as the one indicated by the 
data. 

We also mention that the gravitational wave contribution to the spectrum is not expected to 
be large. For our choice of the parameters of the model (see next section), it is more than two 
orders of magnitude smaller than the scalar contribution during the second stage of inflation. 
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Due to the similarity of the potentials of the first and second stage, we expect that the same is 
true for the first stage of inflation as well. 

Finally, we consider the effect of the subsequent inflationary stages on the spectrum. The 
presence of a third stage, generated by an additional term ^3 (— /^l + As^s^&s) in the superpo- 
tential, could affect our previous discussion. In order to be more specific, let us concentrate 
on the second and third inflationary stages. In the presence of additional superfields $3, $3, 
the superpotential associated with the second inflationary stage has the more general form 
"5*2 (—^2 + X2^2^2 + g'^s^s)- For sufficiently small coupling g', the potential that generates a 
slope along the fiat direction during the second stage has a form analogous to eq. ( |3.11| ), i.e. 



\ 3 



^ Ms 



167r2 



(A3^i + g'fJ'l 



I 2Ai / 2 



. (5.14) 



For simplicity we assume that the value of (T3 does not change substantially during the second 
stage. This is expected to be the case for fj,2 sufficiently larger than /i3. The contribution of the 
fluctuations of the (T3 field to the spectrum generated during the second stage can be estimated 
through the analogous of eq. (|5.8[). These fluctuations do not affect the spectrum if 



/Vc' « 1, (5.15) 

where 

^ M3XI a2i ' 

c=l + /4- (5-16) 

A3 /^3 



The above constraint has implications for the total number of e- foldings N^tot, given by Eq. ( 4.11 ) 



that can be generated during the third stage. For 02 f ^ (72i, csf a^i, we can express N^tot as 

Nstot = N2totf — . (5.17) 

Large values of N^tot require values of / that may be in conflict with the constraint ( |5.15| ). For 
g' = 0, N2tot ~ 2 — 3, cr2i/mpi = 0.025 (the value of this parameter that we use in the next 
section), and allowing for ~ a value of a^i slightly below mpi can lead to the required 
remaining ~ 50 e- foldings. A non-zero value of g' leads to c > 1 and permits even larger values 
of / and, therefore, smaller values of a-^i. 

The discussion of the first stage can be carried out in an analogous way, even though the 
constraints are more complicated. We expect that, similarly to above, there are regions of 
parameter space for which the primordial spectrum of the first stage is not affected by the 
presence of the additional fields. Moreover, it is conceivable that several similar inflationary 
stages generate a total number of e-foldings ~ 50, in which case the resulting bounds are less 
stringent. As we have discussed in the introduction and section 2, we have in mind a picture of 
multiple bursts of inflation, driven by different sectors of the theory, with a sequence of scales 
starting near nipi and continuing down to the scale implied by COBE, or even to lower energy 
scales. In this work we present only a simplifled picture of two of these stages that have direct 
observational consequences. 

Another possibility is that the later stages of inflation are associated with fields that are 
massive during the first stages or located near the minima of their potential, and, therefore, 
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do not affect the spectrum. These fields are displaced from their initial position during the 
reheating and trigger inflation only when the temperature becomes sufficiently low. For such 
models the constraints discussed above do not apply. The scenario of thermal inflation pO| ] is a 
typical example. 

The observational consequences of the third inflationary stage are difficult to estimate. One 
could try to compute the linear mass fluctuation (7{R) for ~ 1 h~^Mpc and compare with the 
abundance of galaxies at large redshifts |44|. However, the spectrum for the relevant region of k 
is strongly affected by the intermediate stage between the second and third inflationary stages, 
for which a computation is impossible. The fact that the scale /is is largely unconstrained by 
previous considerations makes it probable that an agreement with observations is feasible. 



6 Comparison with observations 

In this section we compare the predictions of our model with the experimental data from COBE 
and the APM galaxy survey. We do not consider the gravitational wave contribution to the 
spectrum, as we expect it to be small. We also concentrate on the adiabatic scalar contribution 
given by Eqs. ( |5.10| )-( |5.13 ). We do not consider any early ionization scenarios or HDM. We are 



interested in examining whether the COBE-DMR measurements can be made compatible with 
the data of galaxy surveys employing the smallest number of cosmological parameters. More 
specifically, we consider two cosmological scenaria: 1) A CDM model, with zero cosmological 
constant, present Hubble parameter h = 0.5 {Hq = h 100 km/s/Mpc) and flat geometry 
^matter = 1- We assume that the fraction of the critical density in baryons is = 0.05 and 
in cold dark matter ^cdm = 0.95. The helium abundance is Yhc = 0.24, while we do not 
consider any massive neutrinos. 2) The same model with I^a = 0.5, h = 0.5, = 0.05 and 
^CDM = 0.45. 

In both the above scenaria, we use the following parameters for the model of Section 3: 
Ai = 1, A2 = 0.1, 5 = 7.1 X 10-3, ^i/mpi = 3.9 x 10-^ ^^2/mpl = 8.8 x 10"^, au/mpi = 0.36, 
C2j/?Tipi = 0.025. For the purpose of this section, the parametes Mi^2^ which determine the 
dimensionalities of the representations of the groups Mi^2 to which the superfields $1^2 belong, 
can be absorbed in a redefinition of Ai, A2, g- For this reason we set Mi = M2 = 1. For the 
above choice of parameters the observable part of the first stage of infiation generates Nnot — 5 
e- foldings and a spectrum with a spectral index n ~ 1. For the second stage of infiation, 
Eq. (4.10) predicts N2tot — 2 e-foldings. However, a more accurate numerical integration of 



the field evolution equations gives N2tot — 3. The spectral index relevant for the scales that 
cross outside the horizon in the beginning of the second stage of infiation is given by Eq. ( |5.13 ) 



with (T2{k) ~ o"2i and is n ~ 0.6. The intermediate stage affects the scales k2 < k < ki 
with ki/k2 — 4: A according to Eq. ( |5.3D . The values of the power spectrum for k2 and ki are 
determined by Eqs. ( |5J0| ), (|5.1lD and satisfy P{k2)/P{ki) ~ 11.8. 



The initial spectrum of fiuctuations arising in our model generates matter and radiation 
perturbations, which, after amplification, give rise to the observed large scale structure and the 
anisotropics of the cosmic microwave background. In the following, we compare the theoretical 
predictions of our model with the measured anisotropics of the cosmic microwave background 
and the observed clustering in the galaxy distribution. For the comparison with the CMB data, 
we assume a form of the power spectrum in the range k2 < k < ki that interpolates smoothly 
between the spectrum in the ranges k > ki and k < k2- 
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6.1 Angular power spectrum of CMB anisotropics 

The cosmic microwave background (CMB) radiation, last scattered at the epoch of decouphng, 
follows a blackbody distribution to high accuracy | }45[ ], with a temperature almost independent 
of direction, T = 2.728 it 0.002i('. As it was measured by the DMR experiment on the COBE 
satellite, the CMB radiation has a tiny variation in intensity at fixed frequency, equivalently 
expressed as a variation AT in the temperature: AT/T < 10~^. The 4-year COBE data is 
fitted by a scale-free spectrum, with spectral index n = 1.2 it .3 and a quadrupole normalization 

In Fig. 2 we compare our theoretical predictions for the angular power spectrum of CMB 
anisotropies against the most recent CMB flat-band power measurements (we have used the 
data of Table 1 in Ref. [^]). The coefficients Ce correspond to the expansion of the angular 
correlation function in powers of the Legendre polynomials Pi. For large i, the angular size (in 
radians) of a feature in the sky is related to the order of the multipoles that dominate it by 
£ ~ 1/9. The angular power spectrum of CMB anisotropies is normalized to the four- year data 
from COBE-DMR Eq]. In order to calculate the predictions of our model for Ci we have used 



the code CMBFAST of Seljak and Zaldarriaga |38]. 



In our model, the first stage of inflation leads to a scale-invariant Harrison-Zel'dovich [^] 
(spectral index n = 1) spectrum of perturbations. The intermediate stage and the second stage 
of inflation affect only the highest multipoles in Fig. 2. For this reason, our prediction for the 
first acoustic peak in the case J^a = (short-dashed line) is similar to that of the standard CDM 
(SCDM) model (solid line) with a Harrison-Zel'dovich spectrum at all scales. Only for multipoles 
/ ~ 500-1000 a drop of the values of Ci is observed relative to the standard CDM model. This 



should be contrasted with the double-inflationary model of Refs. [16, 19, 20 1, which predicts a 
first acoustic peak that is too low with respect to the Sachs- Wolfe plateau. This problem can be 
traced to the logarithmic dependence oc ln^/^(i^/A;) of the spectrum for small k in this model. 
Better agreement with the data is obtained for our model with f^A = 0.5 (long-dashed line). 

6.2 Angular catalogues of galaxy positions 

Angular catalogues of galaxy positions provide strong constraints on theories of structure forma- 
tion in the universe. In this subsection we compare the predictions of our model with the power 



spectrum of galaxy clustering derived from the angular APM galaxy survey |30|-|33]. Redshift 
surveys are more noisy than the APM survey at large scales and subject to more uncertainties 
due to the distortion of the pattern of galaxy clustering by the peculiar motions of galaxies. 

As we have discussed in Section 2, the effects of the non-linear evolution of matter fluctuations 
(when the density contrast becomes of order 1) must be removed before a comparison is possible 
with the linear spectrum. This is usually done through A^-body simulations of clustering. Baugh 



and Gaztahaga [33| have estimated the linear spectrum that gives rise to the APM spectrum for a 
spatially flat universe with critical density and zero cosmological constant. Their estimate is well 
fitted by the curve of Eq. (|2.1|) This result is in agreement with the spectrum obtained through 



the application of the formula of Jain et al. |35] for the relation between the linear and non-linear 



spectra. The formula of Peacock and Dodds |36] does not lead to very good agreement with the 



results of the N-hody simulations [ 34 1 . The effects of the non- linear evolution are significant for 
k ^ 0.2/iMpc-^. 

For the comparison of our predictions with the data, we calculate the transfer function T[k) 
to better accuracy than the one provided by Eq. (|2.2D, by employing the code CMBFAST of 



Seljak and Zaldarriaga ||38|]. The linear power spectrum is related to the density perturbation 
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at horizon crossing through [10| 



AHk) = ^^ = 6%ik) fA) r^(fc), (6.1) 



27r2 ^ V^o 

where the dimensionless quantity A'^{k) stands for the mass variance per unit interval in In A: 
{i.e., A'^{k) = dal^,Jdlnk). 

In Fig. 3 we display the power spectrum P{k) in units of (Mpc//i)'^ as a function of k in 
units of /i/Mpc for our model with I^a = and /c2 = 0.06 /iMpc^^. For scales k < k2 the 
spectral index has been taken n = 1, while for scales k > ki the spectrum is tilted with spectral 
index n = 0.6. For scales k2 < k < ki the complicated nature of the intermediate stage makes 
the calculation of the power spectrum very difficult. The prediction of our model for the linear 
spectrum is given by the short-dashed line. The points with error bars correspond to the data 



from the APM angular galaxy catalogue |31, |3^. The fit of Eq. ( |2.lD for the linear spectrum 
derived from these data is given by the solid line. 

We observe good agreement of our predictions with the linear spectrum expected from the 
APM data. This agreement should persist for the scales k2 < k < ki, as no strong new feature 
in the spectrum is expected in this interval (see discussion at the end of Section 5). We point out 
that we have not assumed any biasing in the galaxy distribution with respect to the underlying 
mass fluctuations (bias parameter bg = 1). Allowing for other values of bg provides an additional 
degree of freedom that can improve the agreement of our predictions with observations. 

The power spectrum for our model with = 0.5 is given by the long-dashed line in Fig. 4. 
We compare again with the estimate of Eq. ( |2.1[) for the linear spectrum (solid line). The 
reason is that for a flat geometry of the universe we do not expect a significant modification 



of this estimate. The formula of Peacock and Dodds |36] for the relation between the linear 
and non-linear spectra predicts only a small difference between the spectra for ^matter = 
JIa = and for ^matter = 0.5, Oa = 0.5. We observe a good agreement between our prediction 
and the linear spectrum in the range k > 0.26 hMpc~^. On the other hand, our prediction is 
significantly above the linear spectrum for k < 0.06 hMpc^^. However, the observational data 
in this last range are dominated by large errors. This makes the estimate for the linear spectrum 
less accurate. 

For values of Q\ above 0.5, the discrepancy between our prediction and the APM spectrum 
at large scales is enhanced. Reconciling the prediction with the data can be achieved either 
through the consideration of a bias parameter different from 1 for the galaxy distribution or a 
spectral index smaller than 1 for the first stage of inflation. 

As a final test of our model, we turn to the calculation of the parameter erg, which is a 
common measure of large-scale clustering. It is defined as the variance of the density field 
smoothed over a radius R = 8 h~^Mpc 

1 

al = —j W^{kR)6l{k)T^{k)k^dk. (6.2) 
Hq Jo 

The "top-hat" smoothing function W{kR) is given by 

'sm(kR) cos{kR) 



W{kR) = 3 



{kRY {kRf 



(6.3) 



For n ~ 0.6 the largest contribution to the integration of Eq. ( |6.2[ ) comes from A; ~ = 
0.125 /iMpc~^. This corresponds to the tilted part of the spectrum near ki, for which we have 
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used n = 0.6 in the calculation of the results presented in Figs. 3 and 4. However, as Fig. 1 
indicates, the effective spectral index may be even smaller for k slightly below ki. For this 
reason we have estimated erg allowing for a local variation of n in the interval [0.2, 0.8]. We find 
that fjs is not very sensitive to the value of n. Instead, it is determined by the normalization 
of the spectrum at A; ~ ki, which is fixed by P{k2) / P{ki) ~ 11.8 for our choice of parameters 
of the model. The predicted values of ag are near 0.75, both for Q.\ = and = 0.5. They 
are in good agreement with the values deduced from the abundances of rich clusters of galaxies: 
as ~ 0.6 ± 0.2 for ^matter = 1 , = § , and fTg ~ 0.85 ± 0.3 for ^matter = 0.5, = 0.5 i. 



7 Summary and conclusions 

In this paper we considered the possibility of a multiple-stage inflationary scenario within the 
context of supersymmetric hybrid inflation. This framework has the nice feature that the ob- 
servable part of inflation takes place for field values below the Planck scale. This permits the 
discussion of corrections to the tree-level picture (such as supergravity corrections) in terms of 
controlled expansions in powers of fields in units oimpi. Inflationary models such as the one we 
considered have been discussed in Ref. |22|, and have been shown to survive the supergravity 



corrections in Ref. [23| 



One motivation for the consideration of a scenario with more than one stages of inflation 
arises from the discrepancy between the inflationary energy scale implied by the COBE data 
and the Planck scale at which classical general relativity becomes applicable. This discrepancy 
results in the necessity of an extremely fine-tuned field configuration at the end of the Planck 



era for the onset of hybrid inflation [26|. The problem can be avoided if inflation takes place in 
two stages [^2|. The first stage has a typical scale ~ mp; and occurs naturally. The second has 
a characteristic scale much below mpi and generates the density perturbations that result in 
the cosmic microwave background anisotropy measured by the DMR experiment on the COBE 
satellite. It is the logical next step to contemplate the possibility of a multiple-stage inflationary 
scenario, with a sequence of scales between mpi and close or below the one implied by COBE. 

Another motivation for the consideration of a multiple-stage scenario stems from the need 
to reconcile the power spectrum of density perturbations at large scales (obtained through the 
COBE-DMR measurements for the anisotropies of the cosmic microwave background) with the 
observed power spectrum at small scales (deduced from galaxy surveys). The standard CDM 
model with a flat (spectral index n = 1) initial spectrum of adiabatic perturbations predicts too 
much power at small scales, if one normalizes it to the COBE data at large scales. An interesting 
possibility in order to resolve this discrepancy is to consider power spectra with a break. A strong 
indication for the location of the break is provided by the form of the power spectrum of galaxy 
clustering derived from the APM galaxy survey [|l|, H, H, H. At h ~ 0.05 h Mpc-^ the data 



support a sharp drop of the spectral index of the primordial spectrum, from a value n ~ 1 to an 
average value n ~ 0.6. This break has been linked to a possible phase transition between two of 
the inflationary stages in a multiple inflationary scenario [p^]. 

In this paper we presented a model of supersymmetric hybrid inflation that generates a 
spectrum with a break. As we argued in Section 2, the large deviation of the spectral index 
from 1 for k ^ kb — 0.05 h Mpc~^ results in a small number of e-foldings generated by the stage 
of inflation relevant for this range of scales. A picture of multiple short bursts of inflation, with 
significant variations of the spectral index of the primordial spectrum, seems probable. 

We computed the spectrum of adiabatic density perturbations for the range of scales that are 
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relevant for the COBE data and the APM galaxy survey. This requires the explicit discussion of 
the first two observable stages of inflation, which generate ~ 8 e-foldings. The total number of e- 
foldings must be ~ 60 for the flatness and horizon problems of standard cosmology to be resolved. 
However, the stages responsible for the last ~ 50 e-foldings generate density perturbations at 
scales that are strongly affected by the non-linear evolution. As a result, it is very difficult to 
extract information on their characteristics. 

In our model, the adiabatic density fluctuations during both stages of inflation originate in 
the quantum fluctuations of the same field. The first stage of inflation generates ~ 5 e-foldings 
and a spectrum with index n = 1 for scales k ^ 0.06 h Mpc~^. The second stage generates ~ 3 
e-foldings and a spectrum with index n ~ 0.6 for scales k ^ 0.26 h Mpc^^. The calculation 
of the spectrum for 0.06 h Mpc^^ ^ k ^ 0.26 h Mpc^^ is prohibited by a very complicated 
intermediate stage of normal expansion. However, it probable that no significant new feature 
arises in this range of scales, as we discussed at the end of Section 5. 

We considered two scenaria: 

a) A CDM model, with zero cosmological constant, present Hubble parameter h = 0.5{Hq = 
h 100 km/s/Mpc) and flat geometry - ^matter = 1- We assumed that the fraction of the critical 
density in baryons is O;, = 0.05 and in cold dark matter ^cdm = 0.95. 

b) A scenario with = 0.5, Q;, = 0.05 and ^cdm = 0.45. 

The comparison of our predictions with observations is given in Figs. 2-4. For the first scenario, 
the angular power spectrum of CMB anisotropics (Fig. 2) is very similar to the one with a 
Harrison-Zel'dovich spectrum at all scales. Only for multipoles / ~ 500-1000 a drop of the 
values of Ci is observed relative to the standard CDM model. This should be contrasted with the 
predictions of a model of double chaotic inflation, for which the first Doppler peak is significantly 
suppressed [20|. The second scenario leads to a better agreement of the generated spectrum 
with the data. The power spectrum predicted by the first scenario compares well with the linear 
power spectrum deduced from the APM data (Fig. 3). The agreement is less satisfactory for 
the second scenario (Fig. 4), even though the discrepancy appears in the region of large data 
errors. For both scenaria, our predictions for (Tg are in agreement with the values deduced from 
the abundances of rich clusters of galaxies. 

As a closing remark, we point out that the questions of the presence and the nature of a 
feature in the power spectrum of the galaxy distribution are not settled. Several interpretations 
and explanations of the data exist (e.g., see Refs. |17, 0). Our work provides a possible 
explanation of the observed feature within a multiple-stage inflationary scenario in the context 
of supersymmetric hybrid inflation. 
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Figure 1: The logarithmic slope of: (a) The power spectrum recovered from the APM angular 
galaxy catalogue (points with error bars), (b) The estimated linear spectrum that gives rise to 
the observed spectrum (solid line), (c) The estimated primordial spectrum (dashed line). (After 
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Figure 2: Theoretical predictions for the angular power spectrum of CMB anisotropics in our 
multiple-stage inflation model, against the most recent CMB flat-band power measurements. 
We plot \/t(Jl -\- l)Ct,/{2'K) in units of jiK versus the multipole moment i. The solid line corre- 
sponds to the standard CDM model = 0, n = 1). The short-dashed and long-dashed lines 
correspond to the spectra in our model for f^A = and Q,^ = 0.5 respectively. 
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Figure 3: Power spectrum generated in our model (short-dashed line) for I^a = 0, plotted 



together with the APM power spectrum (data points) and the linear spectrum of Eq. (2A) (solid 
line). The spectrum is flat for k < 0.06 /iMpc^^ and tilted with n = 0.6 for k > 0.26 /iMpc~^. 
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Figure 4: Same as in Fig. 3 for ^\ = 0.5. 
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